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Abstract: 



In this article, we compute the topological expansion of all possible mixed-traces in a 
I hermitian two matrix model. In other words we give a recipe to compute the number 



of discrete surfaces of given genus, carrying an Ising model, and with all possible 
given boundary conditions. The method is recursive, and amounts to recursively 



^ . cutting surfaces along interfaces. The result is best represented in a diagrammatic 

CN ■ way, and is thus rather simple to use. 
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1 Introduction 

X 

; 1.1 Counting surfaces with given boundary conditions 

The problem of boundary conditions is a very important one in statistical mechanics, 
conformal field theory, string theory. . . (see for example [21 [ISl [IH] for recent develop- 
ments). In this article we address the problem of counting configurations of an Ising 
model on a random lattice, with given boundary conditions. This problem can be 
equivalently stated as computing mixed traces expectation values in a 2-matrix model. 

The 2-matrix model was introduced by Kazakov [T^ as the Ising model on a random 
lattice. Its partition function reads: 

dMidM2 ^-NTr[V,(M,)+V,(M,)-M,M,] 
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Figure 1: Example of surface generated by (Tr Mf Tr M|)^, where one has associated 
the blue color to Mi and the red color to M2: it is a cylinder with one boundary of 
length 5 with red condition and one boundary of length 3 with blue condition. 

where Vi and V2 are polynomials, and where the integral is a formal hermitian matrix 
integral (see for example [13] for a definition of formal integrals), i.e. it is computed 
by first expanding the exponential of the non-quadratic part of Vi and V2, and then 
exchanging the sums and integrals. A formal integral is thus a formal series whose 
general terms are moments of gaussian integrals [13] . 

It is well known from Wick's theorem that such a formal integral is a combinatorial 
generating function which enumerates discrete surfaces (also called maps in the com- 
binatorists litterature) whose faces can have 2 possible colors 1 or 2, or let us say + or 
— , or blue or red. 

The moments: 

< Tr M| > (1-2) 

are generating functions for discrete connected surfaces with one boundary of color 1 
and length I (more precisely, surfaces with one marked face of color 1 and of degree 
/, and one marked edge on the boundary, removed from a closed surface). Similarly, 
< Tr M2 > is a generating function which counts surfaces with one boundary of color 2 
and length /. More generally, < Tr Tr M{^ . . Tr M|'" Tr Tr M^" . . . Tr M^""' >^ 
is a generating function which counts connected surfaces with m boundaries of color 1 
and respective lengths h, . . . ,lm, and m' boundaries of color 2 and respective lengths 
/'i,...,/^/ (see fig iH for an example). The subscript <>c in the expectation values 
means "connected part" or "cumulant", it ensures that only connected surfaces appear 
in the Wick expansion. 
More interesting is: 

< TrMjM^' > . (1-3) 

It is a generating function which counts surfaces with only one boundary of length 
/ + 1', with / color 1 sites followed by I' color 2 sites (see fig|2]for an example). 
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Figure 2: Example of surface generated by ( Tr M^M|)^: it is a disc with one boundary 
of length 7 with red condition for 5 adjacent segments followed by two segments with 
blue condition. 

And more generally, 

< Tr M|i mJ M['m} ...> (1-4) 

counts surfaces with one boundary of length Yl^i ~^ K with li sites of color 1 followed 
by l[ sites of color 2 then I2 sites of color 1, etc. 

It is easy to see that one can design such observables for any given boundary 
conditions: any number of boundaries, and any pattern of sites on the boundaries. 

In this article we show how to compute those generating functions for surfaces of 
given topology. 

1.2 Outline and main results 

The paper is organized as follows: 

• in section 2, we summarize briefly some previous knowledge of formal 2- matrix 
model integrals. Namely, we recall how to compute the "disc amplitude", and 
the spectral curve, and from there the result of [H], i.e. how to count surfaces 
with uniform boundary conditions. 

• in section 3, we define appropriate notations for describing arbitrary boundary 
conditions. We recall which cases were already known in the literature. 

• in section 4, we give the formula for computing the generating functions counting 
surfaces of any genus and arbitrary boundary conditions. The formula is best 
represented diagrammatically, and has a very intuitive interpretation. 

• in section 5, we show some examples of applications of our formula, and in par- 
ticular we show how to recover previously known cases. 
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• section 6 is the conclusion. 



• the proof of the main formula of section 4, is written in the appendix, because it 
is rather technical. 

2 Reminder 2-matrix model 

The 2-matrix model has generated a considerable number of works. Here, we use 
the method of loop equations [HI El [20], which is well suited for genus expansion 
computations. 

2.1 The resolvent 

The resolvent is defined as: 

^.W = (t^^) = E^(t^^^!> (2-1) 

^ ^ ' l=Q 

it is a generating function for a disc of color 1 (i.e. discrete surface with only one 
boundary of color 1 and of length /), and x is a complex "fugacity" conjugated to the 
boundary length jf . 

Like any expectation value in a formal matrix model [31 [6] , it admits a topological 
expansion: 

oo 

{x) = J2 ^? i^) ^'"'^ (2-2) 

9=0 

where W^f\x) is the generating function for discrete surfaces of genus g. 

The loop equations which allow to compute W^f^ have been known for a long time 
[20]. More recently, W[^^ was computed for any g in [1 [ni [n]. The result for wf^ 



can be written in terms of an algebraic equation. Let: 

y{x) = Vl{x)-Wf{x). (2-3) 
y{x) is solution of the following algebraic equation [TJ [8]: 

= E{x, y{x)) = {V;{x) - y{xmVM^)) - x) - p(°)(x, y{x)) + 1 (2-4) 

where 

and y must be chosen as the branch of the solution of E{x,y) = which behaves like 
V{{x) for large x. 



^Remark that these resolvents are properly defined when the fugacity x ^ oo 
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2.2 The spectral curve 

In general, correlation functions are multivalued functions of x, and it is better to write 
them as functions on a Riemann surface. 

Therefore, we view x and y as two meromorphic functions living on a compact 
Riemann surface S. 

E{x, y) = ^ 3p E T, / X = x{p) and y = y{p) (2-6) 

Since the equation E{x, y) = has deg V2 solutions in y for a given x, it means that 
for every point p in E, there are deg V2 points in S such that: 

= 0, . . . , ci2 , x{p') = x{p) (2-7) 

where d2 = deg V2, and by convention we assume p^ = p. 

Similarly, if we regard x as a function of y, then the equation E{x, y) = has deg Vi 
solutions for a given y, which means that for every point p in E, there are deg Vi points 
in S such that: 

Vz = 0,...,cii , y{f)=y{p) (2-8) 
where di = deg V/, and by convention we assume p^ = p. 

2.3 Examples 

• If the algebraic curve S build from E{x, y) = has genus zero, it is possible to find 
a rational parametrization [7[ [5], i.e. x{p) and y{p) are rational functions of p: 

X{p) =-fp+ E£o^' /9 Q\ 

yip) = IP + Efcio /^fcP 

where the coefficients a^, /^^ and 7 are determined by y{p) ~p^oo " V^(p) + 

0(p-2) and x{p) ~p^o - ^/vip) + 0{p^). 

In that case the compact Riemann surface S is the Riemann sphere. 

This is the case which counts the Ising model bicolored maps. 

• If the algebraic curve S build from E{x, y) = has genus 1, it is possible to find 
a parametrization with elliptical functions. 

Spectral curves E{x,y) = of genus g > 0, are not generating functions which 
counts maps, but they are still solutions of the loop equations, they have a more 
complicated combinatorical interpretation, and they are very useful for applications to 
string theory for instance. In what follows, we assume that the spectral curve may 
have any genus, and one should keep in mind that only the genus zero case really 
corresponds to the Ising model on random surfaces. 
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3 Definitions 



We assume that the spectral curve E{x,y) = is known, and that x and y are two 
meromorphic functions on the compact Riemann surface E. 

3.1 Notations 

The most general boundary condition for a discrete surface generated by the 2-matrix 
model is made of several boundaries, some of them having color 1, some having color 
2, and some having mixed color boundaries. 
Let us say that we have: 

• m boundaries of color 1, with conjugated parameters x{pi), . . . , x{pm), 

• 11 boundaries of color 2, with conjugated parameters . . . , y{qn), 

• I mixed boundaries such that the i*^ boundary is made of 2ki changes 
of colors. It can be parameterized with 2ki conjugated length parameters 
[x(Pi,i),y(qi,i),x(pi,2),y(qi,2),x(pi,3),y(qi,3), x(pi,k),y(qi,k)]- 

Notice that the pi's and q/s are points on the curve E. 

The generating function for discrete surfaces with that boundary condition is: 



Hki,...,ki;m;n{Si, S2, ■ ■ ■ , Sf, Pi, . . . , Pm', • • • , ?n) 




+Sl,oSm,lSn,o{y{pi) - Vl{x{pi))) + 5ifi5m,oSn,l{^{qi) " "^^2 



(3 



1) 



where 




Tr 



( 



x{pi^i) - Ml y{qi,i) - M2 x{pi^2 



1 1 



) - Ml y{qi,2) - M2 



1 1 



y{.qi,ki) - ^2 



1 




and 



Si — bi,!' qi,l,Pi,2, qi,2,Pi,3, 9i,3, • • • ,Pi,k, qi,k] 



(3-3) 
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is the ordered set of points {pi^i, qi,i}i=i...k up to cyclic permutations, i.e., using a graph- 
ical representation 




(3-4) 



Each p variable stands for a piece of boundary of color 1, whereas each q stands for a 
piece of color 2. 

Each Hk^^,,,^ki;m;n aduiits a topological expansion: 



TT _ nT2-2g-l-rn~n Trig) 

-nki,...,ki;m;n — -'^ ki,...,ki;m;r. 

9=0 



(3-5) 



where Hjf^ ki-m-n generating function for discrete surfaces of genus g with the 

same boundary conditions (indeed, the Euler characteristic of a surface of genus g with 
I + m + n boundaries is x = 2 — 2(7 — / — m — n). 

We represent H^^ j,^.^.^ graphically as a connected surface of genus g, with / cir- 
cular boundaries, and n + m punctures: 



-^£m;n(^L;Pi, . . . gi, • • • , qn) 




(3-6) 



Since the correlation function H^^ ki-m-n '^^es not depend on the order of the traces 
(i.e. one may permute the S'j's), we may choose one of the boundaries (for example 
Si), and draw it on the exterior, and draw the whole surface in the interior of the circle 
Si- Moreover, because of the cyclic invariance of the trace, we may choose a starting 
point on each boundary (for example pi i) by drawing an anticlockwise arrow on the 
boundary from this poinic. 



■^Remember that the boundaries are oriented according to the sequence of points in the traces of 
the correlation functions. 
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Thus, we represent the correlation function H^^ 



(9) 



a surface S, 



ki, ;m;n 



L;m;ra 



;Sl;pi, 



qn) by 



Sl;Pi, • • • ,Pm] Qi, ■ ■ ■ , <in) which is a disc equipped with g handles, 



/ — 1 holes corresponding to the / — 1 remaining non homogenous boundaries, m white 
marked points corresponding to the homogenous boundaries of color 1 and n black 
marked points corresponding to the homogenous boundaries of color 2. Note also that 
every non homogenous boundary is equipped with a sequence of white and black points 
representing the sequence of boundary conditions. 



1,2 4i,2 




(3-7) 



Notice that the other boundaries S2, ■ ■ ■ , Si also have a marked edge Pi^i — > g^^i whose 
orientation is opposite (i.e. clockwise) of that of 5*1. 

3.2 Previously known results 

Some cases are already known in the literature: 

• Planar case: all Hj^^ ^^.q.q (i.e. planar surfaces only) were computed in 

• Non-mixed boundaries: all functions with only non-mixed boundaries, i.e. 
K-L,n were computed in [T2| 1^ [TTl 115] . 

• Only one mixed boundary with k = 1: H^^.^^ was computed in [15j. 

• In particular the sphere with one puncture is the resolvent: 

<iob) = Wf{p) = Vlixip)) - y{p) (3-8) 

• In particular the sphere with one bicolored boundary is [SI [9l [7] : 

r(0) /r^ _ E{x{p),y{q)) 



HlU{P,Q}) 



{x{p) - x{q)){y{q) - y{p)) 



(3-9) 



Below, we compute all the other ones. 
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4 Diagrammatic solution 



Here, we show the recipe to compute recursively any Hs^^...^Si\m;n- The proof (which 
rehes on loop equations, and is explained in the appendix is very technical, whereas 
the solution is rather simple and can be written pictorially. 

4.1 In equations 

In equations, the recursive solution of the loop equations (see the proof in appendix) 
can be written: 



J^gg ^l!o;o(Pl.l.gl,fcl)'^^('') ^ 



{ Y.hY.A\JB={2,...,l} Yla=2 Zl/,J-^fci-a+l,kB;m-|7|;n-|J|({f'l-"'^l,"' • • • ) }, ^B] Pm/i; Qn/j) 

^ a;(pi,c,)-x(r) 
^ Z^a=2 x{pi,a)-x{r) 

-^a^-l!ii-a+l,kL/{i};m;n({^ > • • -Pl.a-l) ^l.a-l}) {Pl,a, Q'l.a, • • • Pl.fei , Q'l.fei } , Sl/{1}; Pm; Qn) 

_|_ V'^' 1 V 

^-^fc!+fe„kL/{i,i};m;n({'S'l(''^),Pi,Q:, • • • , QiM^PiA^ ■ ■ ■ yPi,a-l, (li,a-l}, SL/{l,i}; Pm; Qn) 

+ J2h Eaub={2,...,«} ^£Ia;|/|;|J|('^i(^)' Pi; qj)^i37™-|/|+i;„_|j| (Sb; r, Pm/i; Qn/J) 
+ ELi H^o%ir)Hjf~'^l^.^.JSi{r), S2,..., Si;pi, . . . gi, . . . , g„) 

+^kL7m+i;n(SK(r); r, Pm; Qn)} 



(4-1) 

It looks terrible, but each term can be represented diagrammatically, and it is in fact 
rather simple and intuitive. Let us notice for the moment that this formula involves 
residues (i.e. contour integrals on S) at various points, in particular the ql which are 
defined in eq J2-8l and were we mean j 7^ 0. 

This formula also involves the function h[^^.q which is given in eq J3-9[ 
All the other terms in the RHS of eq ]4-ll are either some Hg^.l^.Js computed recur- 
sively by the same formula, or some Hq^I^.^ which were computed in [121 111 [HI [IS] ■ 
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4.2 Diagrammatic representation 

It is more convenient to represent equation 14-11 diagrammatically: 




where we explain the meaning of those graphs below. 
4.2.1 Cutting surfaces 

Consider a connected surface S with at least one boundary (i.e. / > 1): 




Let Cut (5) be the set of all topologically inequivalent possibilities of cutting the 
surface along a line pi^i Pi^a (we allow the closed line {i, a) = (1, 1)). When we cut 
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along such a line, we can either get a connected or a disconnected surface. The only 
possibility of getting a disconnected surface is if the point pi^a belongs to S'l, i.e. i = 1, 
and if there is no handle going above the cut. 

Here is the algorithm to construct Cut (5): 

• one first has to choose any ending point pj q, on a mixed boundary and draw a 
path going from the left of the starting point to the left of the ending poiniQ: 

— This point can belong to a boundary different from the starting one, i ^ 1: 



There are ki — 1 such possibilities. 
^The orientation is seen from the point of view of an observer living on the upper side of the disc. 




(4-4) 



There are Ylii=2 ^« possibilities. 

It can belong to the same boundary, i = 1, a 7^ 1 : 




(4-5) 
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— It can be the same as the starting point i — l,a — 1: 




(4-6) 



There is only one such possibihty. 

• Once this ending point is chosen, it remains to fix the position of the handles and 
the other boundaries and punctures with respect to this path. The number of 
inequivalent possibilities depends on the respective position of the starting and 
ending points: 

— If the starting and ending points do not belong to the same boundary, the 
surface is not disconnected by the cut, and every choices are equivalent 
since the left and right side of the path belong to the same component of 
the surface: 



There is only one possibility for the boundaries, punctures and handles 
configuration. 

— If the starting and ending points belong to the same boundary, two different 
configurations can occur: either the path does cut the disc into two discon- 
nected parts, i.e. no handle goes above the path. In this case, one has to 
choose for each handle and boundary whether it lies to the left or the right 




(4-7) 
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of the path: 




(4-8) 



There are 2"+"*+^ such configurations; 

Either the path does not separate the disc into two parts, and all the posi- 
tions of handles, punctures and boundaries are equivalent: 




(4-9) 



There is only one such configuration because one can transport the handles, 
punctures and boundaries across the handle above the path. 

We have then built the set Cut(»S) of cut surfaces associated to any surface S. 



4.3 Weights of graphs 

Now, let us associate a weight to each cut surface. We define recursively a weight V 
on the set of graphs: 

Definition 4.1 The weight V of an uncut surface is given by the corresponding corre- 
lation function: 

^ ('5£m;n(SL;Pl, . . . , Pm', Ql, ■ ■ ■ , Qn)) HI^I^^.^{Si,; Pl, . . . , Pm] Ql, ■ ■ ■ , Qn)- (4-10) 

The weight of the disconnected union of two surfaces is the product of their respective 
weights: 

V(S[jS') := V(S) X V{S'). (4-11) 

The weight of a cut surface is obtained from the weight of the surface(s) obtained 
by cutting along the path 7 following the rules: 
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// the starting and ending points 0/7 do not coincide: 



V 



Res 



-p,g^p' {x{p) - x{r)){y{q) - y{r)){x{p') - x{r))H[%Q{r,q) 



-V 



If the starting and ending points coincide: 



Res 



r— >p 



■9''P' {x{p) - x{r)){y{q) - y{r))H[%Q{r, q) 



-V 



(4-12) 

\ ) 
(4-13) 



where the pi 's are the points encircled inside the closed loop. 

With such notations, equation 14-11 can be reinterpreted as: 

Theorem 4.1 The weight of a given surface is equal to the sum of the weights of all 
corresponding cut surfaces: 

V{S)= J2 ^(^) (4-14) 

5Gcut{5) 

I.e. graphically: 

'•i.k, r n 

,i \ - 





+ 






g-h 


Pl.a 


1l,a-l 












+ 




+ 




(4 



Performing this procedure recursively on any correlation functions, one can elim- 
inate the mixed boundaries step by step until there is no mixed boundary left, i.e. 
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until there are only punctures left. The correlation functions with only punctures are 
computed in [T2ll^[TT]. 

5 Examples of applications 

In this section, we show how to use our formula to recover some previously known 
results, in particular the planar case, and surfaces with uniform boundaries. We also 
compute two simple examples: the generating function of discs with four boundary 
operators and the generating function of cylinders with two boundary operators on 
each boundary. 

5.1 Link with former results 
5.1.1 Planar mixed traces 

If one is interested in the planar mixed correlation functions with only one boundary, 
the recursion relation simplifies to: 




One can thus draw the result of the whole recursive procedure as the sum over all 
possible link patterns on the starting disc in such a way they separate all boundary 
variables. This reproduce the decomposition used in [H] to compute the building 
blocks Fk = Cfrf. 

Example: 

The three point mixed correlation function reads: 
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which gives: 



<o(Pi,^2,P.,g.,P3,g3) 



AO) 



f(0) 



(0) 



r(0) 



{^{Pl)-^{r)){y{q3)-yir)){x{p2)-x{r))H[°lg{r,q3){x(p2)-x{r')){y{q3)-y{T')){x{^^^^ 



+ Res 



Res 



<0;o(Pl.g3)<(];o(r,g2)<o'^o(p3,g3)gl°^;o(p2,92)<o^^o(r',gl) 



AO) 



AO) 



{xipi)-x{r)){y{q3)-y{r)){x{p3)-x{r))H[°lg{r,q3){x{r)-x{r')){y{q2)-y{r')){x{p2)-x{r'))Hl"^^^^ 

(5-3) 



(0) 



One can easily show that this coincide with the result of [m] by using explicitly the 
orientation-reversing symmetrj{§ of the correlation function: 



(0) 



(5-4) 



5.1.2 Simple traces topological expansion 

One can remark that all this recursive procedure supposes that the non-mixed correla- 
tion functions are known, since this new diagrammatic representation does not allow to 
compute them. Nevertheless they were computed by a similar procedure in [121 IH E] 
in terms of trivalent graphs and we show that these former rules could be written in a 
graphical representation similar to the one presented in this paper. 

Let us represent VF^^!|^\(p, pi, . . . ,Pk) '■= -f^o'l+i oiP^Pi, ■ ■ ■ iPk) as a disk with k punc- 
tures instead of a sphere with k + 1 punctures (we have drawn the surface generated 
by this function inside the boundary corresponding to p): 





(5-5) 



The recursion relation of jH [TT] 





+ p- 



- -/(h-i) \^ 



(5-6) 



^Combinatoricaly, this means that summing over all oriented surfaces is equivalent to summing 
over all surfaces with the orientation reversed. 
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can then be written: 






• 


• 





+ 




(5-7) 



where once again, the different terms on the RHS are obtained by drawing a basis of 
homologically independent paths on the disk starting and ending on the boundary, and 
the weight of a cutting along this path follows: 

P 




V Res 



dE,{p) 



iviQ) - yi(l))dx{q) 




(5- 



where one sums over all branch points and q is the point conjugated to q (see [H [TT] 
for more details). 



5.2 Four point function on the disc 

The correlation function i^^^HolPi; 1i}P2, (I2) has already been computed in [8l[Tl]. Nev- 
ertheless, this computation used an Ansatz and a symmetry property of the correlation 
function explicitly. Let us recover the same result without using any symmetry con- 
sideration, but using our recursive formula instead. 
The solution of the loop equations reads graphically: 

P, 5i 




(5-9) 



which is translated intc0 
Hi'\p„q„p,,q,)= Res 



h[^^ {Pi, (12)^1"' {p2, ^2) dx{r) H"^' (r, gi 



r(0) 



(0), 



Pi,P2,g1 Hf\r, q2){x{pi) - x{r))iy{q2) - y(r)) <P2) - x{r) ' 

(5-10) 



^For shortening the notations, we write all along this section H^\pi,qi,p2,q2, ■ ■ ■ ,Pk,<lk) 
Htlo{pi,qi,P2,q2, ■ ■ • ,Pk,qk)- 
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Writing 



one gets 



1 



yir) -y{qi) 



+ 



1 



y{(i2) - y{r) {y{q2) - y{qi)){y{q2) - y{r)) y{q2) - y{qi] 
Hi'\puquP2,q2) 



(5-11) 



Hi^\pu q2)H'l'\p2, q2) dx{r 



(0)/ 



+ 

(5 - 12) 
Since 



Res 
Res 



Hf\r,q,) 



g1 Hf\r, q2){x{pi) - x{r)){y{q2) - y{qi)) <P2) - x{r) 

Hf\pi, q2)Hf\p2, q2) dx{r) {y{r) - y{qi))H[°\r, qi) 

'i {x{pi) - x{r)){x{p2) - x{r)){y{q2) - y{qi)) (y(q^) - y{r))H[°\r, q^) 



Hf\p,q){y{q) - y{p)) 



E{x{p),y{q)) 



(5-13) 



x{p) — x{q) 

the integrand of the second term in the RHS is a rational function in x{r) and it is 
easily checked that the integration contour encircles all its poles (this function is regular 
when x{r) oo). Thus this second term vanishes. 

The first term has no pole at r = q-^, thus it involves only simple poles when 
r — > Pi , P2 and we recover the known result [8], [H] : 

Hi'\pu qr)Hf\p2, g^) - Hf\p,, qi)H['\p2, g^) 



Hf\pi,qi,P2,q2) 



(5-14) 



{x{pi) - x{p2)){y{qi) - y{q2)) 

Even if this new derivation of an old result seems more involved technically, it has 
the advantage of being constructive (the derivation of |T3] was based on an ansatz) and 
does not suppose any additional symmetry of the correlation functions (the derivation 
of [H] was based on the fact that H2^\pi, qi,P2, 92) = Hf\pi, q2,P2, qi))- 

5.3 Generating function of cylinders 

The generating function of cylinders with 2 boundary operators on both boundaries is 
obtained by: 





+ 




(5-15) 



which can be translated into 

Hi%o;oi{Pl^9l}AP2,q2}) 



H[^\pi,qi) dx{r) 



Res 



i H[°\r, qi){x{pi) - x{r)){y{qi) - y{r)) 
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X 

(5 - 16) 



where the second term was computed in [T5 



Hi%oiP2,q2;r) = Res . , \ y^^-^ ' ' ^ ' (5-17) 



Hf\p2,q2)H^X{r;r')dx{r) 



'',r,qV {x{p2) - x{r')){y{r') - y{q2)) ' 



and H'^\.Q{r]r') dx{r)dx{r') is the Bergmann kernel. 



6 Conclusion 

In this article we have found a recursive and graphical method to compute correlation 
functions corresponding to every possible boundary condition for the 2-matrix model, 
i.e. bicolored discrete surfaces. 

The result seems to have a nice combinatorial interpretation, as all the possibilities 
of drawing interfaces (between the + and - spins of the Ising model) in all possible 
ways. However, a combinatorial derivation is missing. 

Also, our result can have interpretations in conformal field theories when one goes 
to the so called double-scaling-limit [HI [5] , and should be compared with recent results 
from Liouville theory fi9[ [T6l [2]. In particular, in [2], our formula for planar disc 
amplitudes is interpreted in terms of the interactions of long folded strings and it 
would be interesting to check the non-planar cases as well. 

It would be interesting also to understand how the structure we exhibit in this 
article, and which seems to be related to integrability like in [H], is related to the 
Langlands programm as claimed by [T6] . 
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Appendix A Loop equations 



Here we prove the equation 14- using loop equations. Loop equations is a standard 
and powerful tool in random matrix theory, they are the Ward identities, or Schwinger- 
Dyson equations, they implement the Virasoro or W-algebra constraints, in combina- 
torics they can be viewed as an extension of Tutte's equations, and in fact they just 
consist in integration by parts, or said differentely, the fact that an integral is invariant 
under (an infinitesimal) change of variable. 

For the 2- matrix model, loop equations were first exploited by Staudacher [SU] , and 
then by many authors, and they led to the solution of [12], HI [15] . 



A.l The loop equations 

In order to prove eq J4-H we consider the change of variables 
5Mi :-- ' 



x{pi,i)-Mi y{pi,i)-M2 x{pi^2)-Mi j/(pi,2)-M2 ' ' ' 3;(pi,fc^)-Mi ?y(pi,fcj)-M2 

nUTr^ 1111 



x(p,,i)-^Mi y(vt,i)-M2 x(pi,2)-Mi y{pi,2)-M2 ' ' ' x{pi^ki)-Mi yipi,ki)-M2 
i=l x{pj)-Mi lls=l y{qs)-M2 ■ 

(1-1) 

Writing that the matrix integral is invariant under this change of variable gives the 
loop equation: 

iy{pi,i) - y{qiM) - ^o^x(pi,i)V'/(x(pi,i)))iffci,...,fc,;m;n(5'i, 5*2, ... , 5'^;^!, . . . , p^; gi, . . . , g„) = 

Yl,A\jB={2 1} Ylii J -f^fci,kA;|/|;iJ|('S'i, ^A] Pi; qj)-f^kB;m-|/|+i;n-i j| (Sr; pi,iPm/i; qn/J) 

1 },Sb 

J^'a-i,kA;UI;|J|({Pi,i''?i.iv--pi,Q-i,gi,a-i},SA;pi;qj)-Ji'Q-i,kA;|/h|J|({Pi.Q''?i.i!---Pi,a-i,gi,a-i},SA;pi;qj) 

X 

El sr^ki 
i=2 2^a=l 



x(pi,a)-x(pi,i) 



x{pi,a)-x{pi,l) 

i}, SL/{i,i}; pm; Qn) 
}, SL/{i,i}; Pm; qn; 

-'^kL;m-l;n(SL;PM/{i};qN)— -'^kL ;m- l;n (^L ;Pm/ {i} ;qN) | 



Pl,l-=Pi,, 



x(pi)-x{pi^l) 



■X 



1 V^fcl 

l^a=2 x'(pi,c)-a;(pi,i) 

-f^Q-l,fcl-Q+l,kL/{i};m;n({Pl,l, 91,1; • • • Pl,a-1) '?l,o-l}j {Pl,a-, Ql,a, ■ ■ ■ PlMi '?l,A.l}) ^L/jl}; Pm; Qn) 

-^Cf— l,A;i~a+l,kL/|i} ;m;7i({Pl,cu 
+ ]77-f^kL;m+l;n(SK; Pm; Qn) 

(1-2) 

where Polxf{x) denotes the polynomial part in x of /, i.e. the sum of the positive 
terms in the large x Laurent expansion of f{x). 
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Let us now write its g order in the topological expansion: 

.,ki;m;ni^^^ ^^2) ■ ■ ■ ^ Si , Pi, . . . , Pm, Ql, ■ ■ ■ , Qn) 



'I2h=l -^o!l;o(Pl,l)-^fci,..!fc,;m;n('S'l, 5*2, ... , 5*;;^!, . . . , p^; gi , . . . , g„ 



ri9-h) 



+ Eh T.A\jB={2,...,i} E/,j ^fci,kA;i/i;iJ|('^i' pi, qj)^^kB;m-i/i+i;n-i ji (^b ; pi,iPm/i; Qn/J 

+ J2h ^A\JB={2,...,1} ^a=2 Y1,I,J -^A:i-Q+l,kB;m-|7|;n-|J|({^'l>"' ^1>"' • • • Pl.^i 5 }' Pm/i! 



x{pi,a)-x{pi^l) 



"T Z^i=2 Z^a=l x(p,,(j)-a:(pi,i) 

},SL/{i,i};PM;qN) 



^kf+ki,ki^/l:^^iy,m;ni{Sl^Pi,a, (li,a, Pi,a+1, ■ ■ ■ , Qi^^PiA^ ■ ■ ■ ^Pi,a-1^ (li,a-l}, SL/{l,i}; Pm; QN; 

-'^kL;m-l;n(SL;PM/{i};qN)| ,^ 



Pl.l-=Pi. 



1 



+ V 

^a=2 x(pi^a)~xipi^i) 



x{pi)-x{pi,i) 
X 



}, {pi,a, qi,a, ■ • -Pi.fci, gi,fci}, Sl/{i}; pm; qn) 
}, • • -Pi.fci, Pm; Qn) 



+-^kl;m+l;n(SK; Pm; Qn)- 

Notice that we have used the normalizations i/o;2;o = (Tr Tr )^ + ^^}^-^2 explicitly. 



(1-3) 



A. 2 Solution of the equations 

We can solve this hierarchy of equations by induction in the number of traces in the 
correlations and the genus. Indeed, one can remark that the RHS of Eq. f ll-3p contains 
correlation functions with either less traces (that is to say less arguments) either lower 
genus compare to the correlation function in the LHS. One also knows that the last 
term of the LHS is a polynomial in a;(pi i) of degree di — 1 and one can compute its value 
in the di points q{^^^ for j = 1 . . . c/i independently of h'^^.^-S^-k'-, Pm; Qn)- 

For this purpose, let us study the behavior of the LHS when pi^i q{j^^- If p lies in 
the x-physical sheet and gi to the ^/-physical sheet, the definition of the correlation 
function reads: 

{y{Pl,l) - y{(llM))Hkl..M;m;SSu ^2, • • • , SuPi, ...,Pm;qu...,qn) = 



I -, rn ^ 

m 



i=l 
I 



s=l 
n 

nTr 



yiqs) - M 



(9) 
2' c 

(9) 



s=l 



yiqs)-M2. 
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^ I ^ m n 

+ + Tr ^) n(iV4.,i + Tr -) n Tr . , _ \{ 



1 y») 

(1-4) 
where one notes: 

rr. 1 rr. ^ 1 1 1 1 

rp rp / 

^ ~ Va;(pi,i) - Ml - M2 x{pi,2) - Mi ?/(gi,2) - M2" ' x(pi,fcj - Mi 

(1-5) 

and 



Va;(Pi,i) - Ml y(gi.i) - M2 x(pi,2) - Ml y(gi,2) - M2 ' ' ' x(pi,fcj - Mi y(gi.fcj - Ma/ 

(1-6) 

These terms are monovalued functions as long as the p and q variables stay in their 
respective physical sheets. When qi^k^ belongs to the y-physical sheet in the vincinity 
of oOy, all its images q{ lie in the x-physical sheet in the vincinity of oo^- Thus, this 
expression vanishes for pi 1 — > q[ Hence the Lagrange interpolation formula reads 

r(0) ^ \TTi9) 



{x{piA)-x{r)){y{r)-y{q))Hi%,{r,q,,,J 



(1-7) 

where we have defined: 

Ult..,kr,m;ni^iPl,l)) ■= ^O^x(pi.i) V;'(a;(Pl,l))^fci!...,fc,;m;n('^l. '^2, • • • , Si]Pi, . . . , p^] qi, ■ ■ ■ , qn) ■ 

(1-8) 

Insert this formula into Eq. fll-3p and get: 

rj(g) /o N O -f^l!o;o(Pl,l7'?l,fci) RHS|pi^:=,. 
-"kL-m-n(SL;j3l,---,Pm;gi,---,gn) = ReS ttt 

' ' (a;(pi,i) - x(r))(y(g) - y(r))if|.(j.o(r, gi,fcj 

(1-9) 

where RHS denotes all the terms in the right hand side of Eq. f ll-3p . 

One can simplify some of the terms by changing the integration contour. Indeed, 
consider any term of the form dp-^ (^ x(p/i)-x{p- ) ) RHS of Eq. fll-3p . one can 

compute its contribution to the preceding formula: 

_ ^ j^^g _ ^i°0;o(Pl.l'9l,fcl)/(P'.") 



® This term does not vanish when pi^i qi,ki because of the discontinuity of these functions when 
pi^i changes x- sheets. 
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since one can check that the integrand is a polynomial in x(r). We can now move the 
integration contour in the x basis and we get: 

, -D i^j°o;o(Pi,i,gi,fci)/(Ka) 

dp^^ Res . (1-11) 

- x{r)){y{q) - 2/(r))(x(pi,i) - x{p^,a))H\.6■oir, 

This residue can be evaluated by using one more time Eq. (11-31) and recalling that only 
the terms of the form -ffo;2;o(^)Pj,a) have such poles. This finally gives the result, i.e. 
eqlHl 



-^kl'';m;n(^L5 Pi' • • • ' Pm] Ql, ■ ■ ■ , Qn) — 

{ ELi ^S;o(^)^t.tl;m;n('5i(^)> ^2, ■ ■ ■ , Si;pi, . . . , p^; q^, . . . , g„) 

+ Eh Eaub={2,...,«} E/,j ^fei,kA;|/|;|j|('^i(^)' ^a; Pi; qj)^iB;m_|/|+i;„-|j|(SB; r, pm/i; qn/J) 

+ E/i EylUB={2,...,0 Ea=2 E/,J -^fc^'-a+l.kBim-l/ljn-l J| (fe."' ^1."' • • • , ?l,fci } , Sr; Pm/I; QN/J; 

a;(pi,(i)-a;(r) 
^i=2 Z^o=l x(pi_Q,)-a;(r) 

^-^fci+fe„kL/{i_i};m;n({'5'l('^),Pi,a5 (li,a, Pi,a+l, • • • , qi,ki,Pi,l, ■ ■ ■ ,Pi,a-l, ^i.a-l}? SL/{l,i}; Pm; ^.n) 
^ Z^a=2 x(pi,,)-x(r) ^ 

-^a^-l!li-a+l,kL/{i};m;n({^> • • -Pl.a-b {Pl,a, Qha, ■ • • Pl.fci , ^l.fci } , Sl/{1}; Pm; Qn) 

+^£Tm+i;n(SK(r); r, pm; 



(1-12) 

This recursion equation is a triangular system, thus it allows to compute any H{S) 
recursively. 
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